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The transverse motion inside a Stark decelerator plays a large role in the total efficiency of 
deceleration. We differentiate between two separate regimes of molecule loss during the slowing 
process. The first mechanism involves distributed loss due to coupling of transverse and longitudinal 
motion, while the second is a result of the rapid decrease of the molecular velocity within the final 
few stages. In this work, we describe these effects and present means for overcoming them. Solutions 
based on modified switching time sequences with the existing decelerator geometry lead to a large 
gain of stable molecules in the intermediate velocity regime, but fail to address the loss at very low 
final velocities. We propose a new decelerator design, the quadrupole-guiding decelerator, which 
eliminates distributed loss due to transverse/longitudinal couplings throughout the slowing process 
and also exhibits gain over normal deceleration to the lowest velocities. 

PACS numbers: 33.55.Be, 39.10.+j, 39.90.+d 



I. INTRODUCTION 

Recent development of cold polar-molecule sources 
promises to reveal many interesting, and hitherto unex- 
plored, molecular interaction dynamics. The permanent 
electric dipole moment possessed by polar molecules pro- 
vides a new type of interaction in the ultracold environ- 
ment. This electric dipole-dipole interaction (and control 
over it) should give rise to unique physics and chemistry 
including novel cold-collision dynamics [1, 2] and quan- 
tum information processing [3]. 

To date, cold polar- molecule samples have been 
produced most successfully via three different mecha- 
nisms: buffer gas cooling [4, 10]; photo- and magneto- 
association [5-8]; and Stark deceleration [9]. Buffer gas 
cooling achieves temperatures below 1 K through ther- 
malization of molecules with a He buffer. This technique 
produces relatively large densities (10^ cm~^) of polar 
ground- state molecules. However, cooling below ~ 100 
mK has not yet been achieved because the He buffer 
gas has not been removed quickly enough for evapora- 
tive cooling [10]. Photoassociation achieves the lowest 
molecular temperatures of these techniques (~100 /iK), 
but is limited to molecules whose atomic constituents 
are amenable to laser- cooling. Furthermore, molecules 
in their ground vibrational state are not readily pro- 
duced, yielding species with relatively small effective elec- 
tric dipoles, although this problem can be overcome with 
more sophisticated laser control techniques [6, 11]. Stark 
deceleration exists as an alternative to these methods 
as the technique employs well-characterized supersonic 
beam methods [12] to produce large densities of ground 
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state polar molecules (~10^ cm~^), albeit at high packet 
velocities. One limitation of this technique for trapping 
of decelerated molecules is an observed drastic loss of 
slowed molecules at very low mean velocities in both our 
own group's work and the Berlin group of G. Meijer [13]. 
We address this problem in this Article. 

As the leading method for producing cold samples of 
chemically interesting polar molecules. Stark deceleration 
has generated cold samples of CO [9], ND3 [14], OH [IS- 
IS], YbF [19], H2CO [2], NH [20], and SO2 [21], lead- 
ing to the trapping of both ND3 [14] and OH [18, 22]. 
Given the importance of Stark deceleration to the study 
of cold molecules, it is crucial that the technique be re- 
fined to achieve maximum deceleration efficiency. In this 
work, we provide a detailed description of processes that 
limit the efficiency of current decelerators and propose 
methods for overcoming them. We propose possible solu- 
tions to the parametric transverse/longitudinal coupling 
loss originally highlighted in Ref. [23], as well as elu- 
cidate a new loss mechanism unique to producing the 
slowest molecules. We restrict the described simulations 
and theory to Stark decelerated, ground-state OH radi- 
cals, as the supporting experimental data was taken with 
this molecular species. However, the loss mechanisms 
described herein are not specific to OH, and represent 
a general limitation of current Stark decelerators. This 
reduces the efficiency of Stark decelerators for trapping 
cold polar molecules. 

This Article is organized as follows. Section II de- 
scribes the mechanisms responsible for molecular loss 
at low velocities. Sections HI, IV, and V present 
methods of producing molecular packets at intermedi- 
ate velocities (>100 m/s) without the distributed trans- 
verse/longitudinal coupling losses of Ref. [23]. However, 
these methods exacerbate the problem of low-velocity 
loss. Therefore, we propose a new decelerator design in 



2 



Section VI that exhibits gain over conventional Stark de- 
celeration to velocities as low as 14 m/s. 



II. LOSS AT LOW VELOCITIES 

In much of the previous work on Stark deceleration, it 
is assumed that all motion parallel and transverse to the 
decelerator axis is stable up to some maximum excur- 
sion velocity and position from the beam center [17, 24], 
enabling the derivation of an analytical solution to pre- 
dict stable-molecule phase-space area. However, there 
are several important instances where the assumptions 
of this model become invalid. In the case of very slow 
molecules (<50 m/s) [31], we observe two distinct phe- 
nomena leading to reduced decelerator efficiency at the 
final deceleration stages: transverse overfocusing and lon- 
gitudinal reflection. Transverse overfocusing occurs when 
the decelerated molecules' speed becomes so low that the 
decelerating electrodes focus the molecules too tightly 
(transversely) and they either make contact with the elec- 
trodes or are strongly dispersed upon exiting the decel- 
erator. Longer decelerators tend to exacerbate this effect 
due to the fact that molecules can travel at low speeds for 
many stages. Nonetheless, there are several motivating 
factors for constructing a longer decelerator. First, longer 
decelerators allow less energy per stage to be removed 
and consequently lead to larger longitudinal phase-space 
acceptance. Second, a longer decelerator may allow de- 
celeration of molecules possessing an unfavorable Stark 
shift to mass ratio. We will discuss critical issues for use 
of such long decelerators for slow molecule production. 

A second low- velocity effect, which we have denoted 
"longitudinal reflection," is a direct result of the spatial 
inhomogeneity of the electric field at the final decelera- 
tion stage. As highlighted in the context of transverse 
guidance in Ref. [23] , the longitudinal potential is largest 
for those molecules passing-in the transverse dimension- 
nearest to the decelerator rods. However, the decelera- 
tor switching sequence is generally only synchronous to 
a molecule on-axis traveling at the mean speed of the 
packet. As a result, when the mean longitudinal kinetic 
energy of the slowed packet becomes comparable to the 
Stark potential barrier at the last stage, molecules off- 
axis can be stopped or reflected, resulting in a spatial 
filtering effect. Furthermore, the longitudinal velocity 
spread of the molecular packet at the final stage, if larger 
than the final mean velocity, can lead to reflection of the 
slowest portion of the packet. It is important that the 
phenomena of overfocusing and longitudinal reflection be 
addressed since molecule traps fed by Stark decelerators 
require slow packets for efficient loading. 

To illustrate these low- velocity effects, the number of 
phase-stable molecules predicted by three-dimensional 
(3D) Monte Carlo simulation is shown in Fig. 1(a) as a 
function of stage number for increasing phase angle [32] . 
We define the deceleration phase angle, (pQ, exactly as in 
previous publications, where (j)o = {z/L)lSO° [17]. The 
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FIG. 1: (color online) (a) Simulations of the phase- stable 
molecule number as a function of stage number in our 142- 
stage decelerator. Note the dramatic decrease in number in 
the last several stages for = 30.43°. This decrease is due 
to transverse overfocusing and longitudinal reflection of these 
slow (14 m/s) molecules, (b) Simulated transverse (trace 1) 
and longitudinal (trace 2) fractional loss rates as a function of 
time within the final three stages at (f)o — 30.43°. The vertical 
dashed lines denote the times of the given stage switches. 



length of one slowing stage is given by L (5.5 mm for 
our machine), while the molecule position between suc- 
cessive stages is denoted as z. We define z = to be 
exactly between two slowing stages, therefore, <po = 0° 
(bunching) yields no net deceleration. Phase angles 
satisfying 0° < < 90° lead to deceleration of the 
molecular packet, while the maximum energy is removed 
for = 90°. The 3D simulation results displayed in 
Fig. 1(a) include both longitudinal and transverse effects. 
The molecules have an initial velocity (viniuai) of 380 
m/s, corresponding to the mean velocity of a molecular 
pulse created via supersonic expansion in Xe. All simula- 
tions and experimental data hereafter possess this VinUiai 
unless otherwise noted. 

As expected, a higher phase angle leads to a smaller 
number of decelerated molecules. However, there is a 
sharp loss of molecules in the last several deceleration 
stages for the highest phase of (j)o = 30.43°. This value 
of 00 produces a packet possessing a final velocity {v final) 
of 14 m/s. This loss is attributed to transverse overfo- 
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cusing and longitudinal reflection. These distinct effects 
are illustrated in Fig. 1(b), which displays the transverse 
(trace 1) and longitudinal (trace 2) fractional loss rate 
of molecules traversing the flnal three slowing stages at 
(/>o = 30.43°, 

final — 14 m/s. The switching time for 
each stage is denoted by a vertical dashed line, which 
is labeled by the corresponding stage number. Longi- 
tudinal reflection of molecules is clearly the dominant 
loss mechanism for the lowest flnal velocity shown in 
Fig. 1(a). Nonetheless, there also exists a non- negligible 
rise in transverse losses at the flnal stage. That is, be- 
cause the molecular beam is moving very slowly in the 
last few deceleration stages, the transverse guiding flelds 
of the decelerator electrodes have a greater focusing effect 
on the molecules (see Eq. 3 of Ref. [16]) and focus the 
molecules so tightly that they collide with a deceleration 
stage and are lost. In the case of our decelerator, this 
leads to loss of 20% of the decelerated molecule num- 
ber between = 30° (50 m/s) and 0o = 30.43° (14 
m/s). Such a dramatic loss is not predicted by analyt- 
ical theory, as the stable phase-space area decreases by 
only <1% over this range of 00 . This number is calcu- 
lated directly after the decelerator is switched-off and is 
thus an upper bound, since experiments employing these 
cold molecules require them to travel out of the deceler- 
ator where transverse spread can lead to dramatic loss of 
molecule number. 

Experimental evidence of this sudden decrease in 
molecule number at very low velocities is given in Fig. 2, 
which displays time-of-flight (ToF) spectra along with 
corresponding Monte Carlo simulation results at vari- 
ous phase angles. The decelerated OH molecules are 
in the weak-fleld seeking \F = 2, mi? = ±2,—) state. 
The flrst two quantum numbers of the state denote its 
hyperflne components, while the third number indicates 
the parity of the state in the absence of electric flelds. 
Figures 2(b) and (c) illustrate progressive time delays 
of the slowed molecular packets from the background of 
non-synchronous molecules. Figures 2(d) and (e) dis- 
play only the decelerated packets, and the vertical axes 
of these plots are magnifled by a factor of 100. The dual 
peaks correspond to molecular packets loaded into two 
successive stages, and the later peaks represent the in- 
tended flnal velocities. The total detected and simulated 
molecule numbers are plotted in Fig. 2(f) as a function 
of flnal velocity, along with the theoretically expected 
decelerator efficiency (dashed line) [17]. Note that the 
sudden population decrease in both simulation and ex- 
perimental results is not reflected in the one-dimensional 
theoretical model, which does not account for the trans- 
verse dynamics or fleld inhomogeneities that cause such 
behavior. This effect is detrimental to the production of 
dense samples of cold molecules. 

To remove the overfocusing effect at low velocities, the 
transverse focusing of the last several decelerator stages 
needs to be reduced. Different types of transverse focus- 
ing elements may be inserted into the deceleration beam 
line to compensate for this phenomenon. This idea is dis- 



cussed in detail in Sections IV- VI. We note that the pro- 
posed solutions, while successful in addressing the detri- 
mental longitudinal/transverse coupling effects, do not 
mitigate the problem of longitudinal reflection at low ve- 
locities. 



III. DISTRIBUTED LOSS 

Coupling between transverse and longitudinal motion 
throughout the deceleration sequence invalidates the as- 
sumptions made in Ref. [17, 24]. The fact that the trans- 
verse guidance of the molecular beam comes from the 
same electrodes that provide the deceleration means that 
the longitudinal and transverse motions are necessarily 
coupled. While this phenomenon is well understood in 
the fleld of accelerator physics [25], it was flrst pointed 
out in the context of Stark deceleration in Ref. [23] . The 
result can be seen in the left column of Fig. 3, where the 
longitudinal phase space of OH packets is shown versus 
increasing phase angle. In Fig. 3, the dark lines repre- 
sent the separatrix, partitioning stable deceleration phase 
space from that of unstable motion as calculated from Eq. 
2 in Ref. [17]. Each dot represents the position in phase 
space of a simulated molecule. In the absence of cou- 
pling between the longitudinal and transverse motions, 
one would expect the entire area inside the separatrix to 
be occupied. Therefore, the structure in these graphs is 
evidence of the importance of the transverse motion. 

In the left column of Fig. 3, the coupling of longi- 
tudinal and transverse motions is responsible for two 
effects [33]: First, in the center of the stable area at 
00 = 0°-near the synchronous molecule-the density of 
stable molecules is less than in the surrounding area. 
This is because molecules that oscillate very near the 
synchronous molecule experience little transverse guid- 
ing. This effect is not dramatic and is only discernable for 
an exceedingly large number of stages. Furthermore, this 
effect is even less important for the increased phase angles 
typically used for deceleration, since for these switching 
sequences the synchronous molecule experiences more of 
the transverse guiding forces than it does during bunch- 
ing. The second effect, which is much more evident, is 
the absence of molecules at intermediate distances from 
the synchronous molecule as shown in the left column 
of Fig 3. This so-called 'halo' is due to parametric am- 
pliflcation of the transverse motion and is similar to the 
effects seen in cold molecule storage rings [24] as well as 
charged particle accelerators [26]. Essentially, the longi- 
tudinal oscillation frequency of a molecule in this region 
is matched to the transverse oscillation frequency, lead- 
ing to ampliflcation of the transverse and longitudinal 
motion and consequent loss. 

There is a compromise between decreasing longitudinal 
phase-space area and increasing transverse guidance for 
increasing 00 . To demonstrate this effect, we decelerate 
molecules to a flxed v final and vary the phase angle used 
to reach this velocity. This is done either by changing 
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FIG. 2: (color online) (a-e) Experimental ToF profiles (upper 
curves) with Monte Carlo simulation results (lower curves) for 
slowing at various phase angles for ViniUai = 406 m/s. The 
vertical axes of panels (d) and (e) have been scaled up by a 
factor of 100. Note the good correspondence of simulation 
and data to velocities as low as 30 m/s (^o = 34.53°). (f) Ex- 
perimental (dots) and simulated (open circles) total molecule 
number as a function of final velocity for the data in panels 
(a-e). The dashed curve is the expected decelerator efficiency 
calculated from the one-dimensional (ID) theoretical model 
of Ref. [17]. 



the voltage applied to the decelerator rods or by modify- 
ing the effective length of the decelerator itself for each 
^0 of interest. The experimental data shown in Fig. 4 is 
the result of varying the voltage applied to our deceler- 
ator rods (squares) and the effective decelerator length 
(circles). Both lowering the decelerator voltage and us- 
ing shorter lengths of the decelerator for slowing requires 
increasing <po to observe the same v final of 50 m/s. We 
are able to effectively shorten the decelerator by initially 
bunching the packet for a given number of stages be- 
fore beginning a slowing sequence. Note that we use 
5* = 3 bunching to remove any transverse/longitudinal 
couplings during these first stages, then switch back to 
= 1 slowing for the remainder of the decelerator. The 
parameter S signifies the mode of decelerator operation 
and is previously defined in Ref. [23] as well as in Sec- 
tion IV. The phase stable region of 5* = 1 for (j)o > 40° 



FIG. 3: Monte Carlo simulation results for the longitudinal 
phase space of decelerated molecules. The left column shows 
0^ = 0°, 13.33°, and 26.67° for S = 1 slowing, while the right 
column shows (/)o = 0°, 40°, and 80° for S = 3 deceleration. 
The factor of three between 5* = 1 and 5* = 3 phase angles 
ensures that molecules have roughly the same final veloci- 
ties. The observed velocity difference at high phase angles 
occurs because the 142 stages of our slower is not a multiple 
of three. Note that, although the S = 3 phase plot is more 
densely populated than that of 5* = 1 at = 0°, its phase- 
space acceptance decreases dramatically relative to 5* = 1 at 
the lowest velocities. All plots are generated using an identi- 
cal initial number of molecules, and therefore the density of 
points is meaningful for comparison. For definition of 5* = 1 
and 5* = 3, see Ref. [23] and Section IV. 



is entirely contained within that of = 3 at = 0°, 
therefore no artifacts from initial velocity filtering are 
present in this data. We observe that, contrary to the 
predictions of the one-dimensional theory [17], a higher 
phase angle can lead to greater decelerator efficiency up 
to some maximum ^o- This is a direct consequence of 
distributed transverse/longitudinal couplings illustrated 
in Fig. 3. At even larger phase angles, the longitudinal 
phase-space acceptance becomes a limiting factor. The 
labels next to each data point correspond to either the 
voltage applied to the decelerator rods (squares) or the 
number of utilized S = 1 slowing stages (circles). Fig- 
ure 4 further illustrates that the transverse/longitudinal 
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FIG. 4: (color online) Experimental results from changing 
the voltage on decelerator rods (squares) and decreasing the 
effective decelerator length (circles). Effective slower length 
is modified by initially operating the decelerator at 00 = 0°, 
5* = 3, then slowing with 5* = 1 to v final — 50 m/s for the 
number of stages labeled. Both curves illustrate that trans- 
verse/longitudinal couplings are strongly dependent on phase 
angle, and have a marked effect on decelerator efficiency. 

couplings outlined by Ref. [23] reduce decelerator effi- 
ciency, and are highly dependent on phase angle. 

The coupling between longitudinal and transverse mo- 
tion is detrimental to efficient operation of a Stark 
decelerator, reducing the total number of decelerated 
molecules. This effect will be even worse for deceler- 
ating molecules with an unfavorable Stark shift-to-mass 
ratio. Fortunately, the transverse and longitudinal mo- 
tions can be decoupled by introducing a transverse fo- 
cusing element to the deceleration beam line that over- 
whelms the transverse focusing provided by the deceler- 
ation electrodes in analogy to the focusing magnets of 
charged-beam machines. This technique also has the ad- 
vantage of providing a larger stable region in the trans- 
verse phase space which further enhances the decelerated 
molecule number. The remainder of this manuscript dis- 
cusses methods of implementing a transverse focusing ele- 
ment to decouple the longitudinal and transverse motion. 
Sections IV and V describe methods that can be imple- 
mented by modifying the timing sequences of present de- 
celerators with moderate success at intermediate v final- 
Section VI presents an improved design for a Stark de- 
celerator that solves this problem of distributed longitu- 
dinal and transverse loss and also reduces the previously 
described overfocusing losses at the final stage. 



IV. DECELERATOR OVERTONES 

The simplest method for introducing a transverse fo- 
cusing element to the decelerator beam line is to let 
the molecules fly through an energized deceleration stage 
without removing the field. In this manner, molecules ex- 



perience the transverse focusing of the entire stage with- 
out their longitudinal motion affected. Traditional longi- 
tudinal phase stability requires the switching of the fields 
to occur on an upward slope of the molecular potential 
energy, i.e., faster molecules are slowed more while slower 
molecules are slowed less than the synchronous molecule. 
Hence, it is necessary to de-sample the bunching switch- 
ing rate by an odd factor (3,5,7...): the so-called deceler- 
ator overtones [27]. For convenience, we define the quan- 
tity S = Vo/v Switch ^ where Vq is the synchronous molecule 
velocity and the switching speed vswUch is given as the 
stage spacing L divided by the switching time-interval. 
Reference [27] considered only the bunching case. In this 
work, we generalize to the case of actual deceleration. 
However, the above definition of S is still valid. That is, 
S is constant despite the fact that both Vq and vswUch 
vary when > 0°. With this definition we see that tra- 
ditional deceleration can be described by 5* = 1, while 
the method of de-sampling the switch rate by a factor 3 
is described by 5 = 3. These two methods of deceleration 
can be seen in Figs. 5 (b) and (c), where their respective 
switching schemes are shown for = 0°. By switching 
at one-third the rate, the molecule packet flies through a 
deceleration stage that is energized and experiences en- 
hanced transverse guiding. 

Longitudinal phase space simulations of = 3 slowing 
at various phase angles are shown in the right column of 
Fig. 3. No structure is present in these plots. Also, the 
region of longitudinal phase stability for S = 3-even at 
= 0°-is reduced compared to = 1. This is because 
the maximum stable velocity, as calculated from Eqs. 2 
and 6 of Ref. [17], depends on the spacing between decel- 
eration stages as and thus, the separatrix velocity 
bound is reduced by a factor of a/3 [34] . Nonetheless, the 
absence of coupling to the transverse motion leads to a 
larger number of molecules for the = 0° case shown 
in the uppermost panel of Fig. 3. In a given decelera- 
tor, S = 3 slowing requires a factor of three higher phase 
angle than S = 1 to reach the same final velocity. As 
illustrated in Fig. 3, this fact severely limits the practi- 
cality of 5* = 3 as a deceleration scheme, as it implies a 
dramatic reduction in velocity acceptance at the highest 
phase angles. 

Figure 6 shows ToF data (upper curves) taken in our 
142-stage decelerator for = 3 with increasing phase 
angle along with the results of 3D Monte Carlo simula- 
tions (lower curves). Qualitatively, the ToF signal is no 
different than for 6* = 1. 

Shown in Fig. 7(a) are experimental decelerated molec- 
ular packets for 5 = 3 (upper) and 5* = 1 (lower) ver- 
sus increasing phase angle. For each successive packet, 
the S = 3 phase angle increases by 10° in the range 
0o = 0°-60°, while the S = 1 phase angle increases by 
10° /3. In this manner, the packets are decelerated to 
roughly the same velocity. There is a slight difference at 
the highest phase angles shown because the total number 
of stages in our decelerator (142) is not an exact multiple 
of 3. In Fig. 7(b) the de-convolved total molecule num- 
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FIG. 6: (color online) Experimental ToF data (upper curves) 
of OH molecules decelerated using S = 3 at varying phase 
angles. The results of 3D Monte Carlo simulations are also 
shown (lower curves). 



FIG. 5: (color online) Deceleration schemes, (a) Potential 
energy shift of polar molecules in the Stark decelerator. The 
dotted (blue) curves show the potential energy shift when 
the horizontal (circular cross section) electrodes are energized, 
while the dashed (red) curves show the potential energy shift 
when the vertical (elongated cross section) electrodes are en- 
ergized. Deceleration proceeds by switching between the two 
sets of energized electrodes. In panels (a)-(c), the thick black 
line indicates the potential experienced by the molecules. The 
empty circles indicate a switching event, (b) Traditional (S = 
1) operation at = 0°. For phase stability, the switching al- 
ways occurs when the molecules are on an upward slope, and 
as such the molecules are never between a pair of energized 
electrodes. Thus, the maximum transverse guiding is never 
realized, (c) First overtone operation (S — 3) at (/)o = 0°. By 
switching at one-third of the S — \ rate, the molecules are 
allowed to fly directly between an energized electrode pair, 
and thus, experience enhanced transverse guiding, (d) Opti- 
mized first overtone operation (S — 3+) at = 0°: Initially, 
the packet rises the Stark potential created by one set of elec- 
trodes. When the molecules reach the apex of this potential, 
the alternate set of electrodes is energized in addition. In this 
way, the molecules experience one more stage of maximum 
transverse guiding for each slowing stage. Note that, to min- 
imize the un-bunching effect, the grounded-set of electrodes 
is switched on when the molecules are directly between the 
energized electrodes. 



ber for each of these packets is plotted versus final speed. 
While the 5* = 3 method dominates over 5 = 1 for small 
phase angles, its effectiveness decreases as the deceler- 
ation becomes more aggressive-by 224 m/s, the S — \ 
molecule number is already larger than that of 5* = 3. 
This behavior is expected since the phase angle used to 



decelerate to 224 m/s is 60°, while the required S = \ 
phase angle is only 20°. Although the = 3 longitudinal 
phase bucket does not exhibit structure, it is so much 
smaller in enclosed area than the 5 = 1 that its total 
molecule number is smaller. The simulation results of 
Fig. 7(c) and the theory results of Fig. 7(d) support this 
description. Figure 7(c) shows that, even at 0o = 80° , 
the transverse loss rate per stage is at all times greater 
for 5 = 1 than 5 = 3. However, the calculated longi- 
tudinal phase-space acceptance of Fig. 7(d) mirrors the 
behavior observed experimentally in Fig. 7(b) when the 
initial points are scaled to the experimental ratio of 2.75. 
This scaling accounts for the aforementioned 'halo' in the 
slowed 5 = 1 packet, which persists relatively unchanged 
over the range of 5 = 1 phase angles used (0^ = 0° — 20°). 
The fact that the theory curves of Fig. 7(d) cross at a 
higher velocity than the data of Fig. 7(b) suggests there 
is increased transverse guiding of 5 = 3 slowing at high 
phase angles. Nevertheless, even with 142 stages of de- 
celeration, 5 = 3 is unfavorable for velocities below 224 
m/s due to reduced stable phase-space area. 

At this point, one would expect that operation at lower 
phase angles would permit realization of the gain pro- 
duced by the 5 = 3 method. This can be accomplished 
by naively using a longer decelerator. A simulation of 
this kind is presented in Fig. 8, which plots the number 
of molecules present after deceleration at = 20° versus 
final velocity [35] . The number next to a data point rep- 
resents the number of deceleration stages used. Initially 
the decelerated molecule number is relatively flat versus 
final velocity. However, after about 500 stages (180 m/s), 
the number of decelerated molecules begins to decrease 
and dramatically falls off after 550 stages (150 m/s). Very 
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FIG. 7: (color online) Comparison of deceleration using S = 3 
versus S = 1. (a) Experimental ToF data of decelerated OH 
packets with S = 3 (top) and S = 1 (bottom). Note the 
factor of three between S = 1 and 5* = 3 phase angles, (b) 
De-convolved, integrated molecule number for S = 3 (trace 
2) and S = 1 (trace 1) for the packets shown in panel (a), 
(c) Simulated transverse loss rate per stage for 5* = 1 (trace 
1) and S = 3 (trace 2) deceleration. As expected, S = 1 
results in larger transverse loss rates throughout, (d) Calcu- 
lated stable longitudinal phase-space area for 5* = 1 (trace 
1) and S = 3 (trace 2), with initial points scaled to the ex- 
perimental ratio of 2.75. The above panels highlight that the 
observed shortcoming of 5* = 3 deceleration is entirely due to 
loss of longitudinal velocity acceptance at high phase angles. 




Final Velocily [m/s] 



FIG. 8: Monte Carlo results for decelerated molecule number 
using S = 3 and 0o = 20° versus final velocity. The number 
next to each data point is the number of stages used. Because 
of transverse overfocusing and longitudinal velocity filtering, 
essentially no molecules survive below 100 m/s. 



few molecules survive below 100 m/s. This is because for 
5* = 3 the decelerated molecules must fly through an 
entire stage while experiencing a guiding force in only 
one dimension [see Fig. 5(c)]. Once the molecules are 
at slower speeds they can spread out in one transverse 
dimension or be over-focused in the other and collide 
with the rods. As the mean kinetic energy of the slowed 
packet becomes comparable to the full potential height, 
the packet can be nearly stopped as it traverses the inter- 
mediate charged stage. This has two consequences: (1) 
longer transit time leading to more intense transverse fo- 
cusing; and (2) velocity filtering of the low-speed packet 
as slower molecules are longitudinally reflected from this 
potential. 

The transverse loss responsible for the extreme drop 
in molecule number for = 3 deceleration also occurs 
in traditional deceleration, but to a lesser degree. Be- 
cause of this decrease in molecule number at low speeds, 
the usefulness of slowing with S = 3 is generally lim- 
ited to experiments that do not require the lowest veloci- 
ties, such as microwave spectroscopy and collision exper- 
iments [28-30]. We note also that our simulations pre- 
dict no low-velocity gain when using slowing sequences 
containing combinations of deceleration at /S = 1 and 
bunching at S = 3. 



V. MODIFIED DECELERATOR OVERTONES 

A natural extension of the above overtone decelera- 
tion is the use of what we have termed a "modified de- 
celerator overtone," denoted by an additional plus sign, 
i.e., S = 3+. Deceleration in this manner is shown in 
Fig. 5(d) for = 0°. In this method, deceleration pro- 
ceeds similarly to conventional S = 3. However, S = 3+ 
sequences yield confinement of the packet in both trans- 
verse dimensions. This is achieved by charging all slower 
rods for the period in which the synchronous molecule is 
between switching stages. In order to minimally disrupt 
the longitudinal dynamics of the synchronous molecule, 
the second set of slower rods is charged only when the 
molecule is at the peak of the longitudinal potential from 
the first rod set. The packet then traverses two charged 
rod pairs before reaching the next switching stage, at 
which point the rod pair that was originally charged is 
grounded. While this does create a slight anti-bunching 
effect, i.e., molecules in front of the synchronous molecule 
gain a small amount of energy, it provides an extra stage 
of transverse guidance in comparison to S = 3. Exper- 
imental results of this method of slowing are shown in 
Fig. 9 for comparison to deceleration using both S = 1 
and = 3. Fig. 9(a) displays a unique consequence of 
the S = 3+ switching sequence, where operation at 0^ 
= 0° leads to deceleration. Figure 9(b) shows that op- 
eration with = 3+ provides slightly more molecules 
than 6* = 3, but the loss of molecules due to decreased 
longitudinal velocity acceptance remains a problem. The 
increase in molecule number for S = 3+ over S = 3 is 
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FIG. 9: (color online) (a) Experimental ToF data of deceler- 
ated OH packets produced using the S = 3+ modified over- 
tone. Also shown for comparison is the experimental bunch- 
ing packet for operation at 5* = 1. (b) The de-convolved, 
integrated molecule number calculated from 5* = 3+ (open 
squares), S = 3 (open circles), and S = 1 (dots) data. 



due to the extra stage of transverse guidance, which for 
the higher- velocity packets we measured leads to a larger 
transverse acceptance. 

To determine whether the extra stage of transverse 
guidance would counter the overfocusing effects [36], we 
perform simulations of 5* = 3+ deceleration at (j)o = 
20° for a varying number of deceleration stages. The 
results of these simulations, shown in Fig. 10, are simi- 
lar to the results for 5 = 3. Namely, as the decelerator 
length is increased and the molecules' speed is reduced, 
there is a marked molecule number loss for velocities be- 
low 200 m/s. In our simulations, we could not observe 
any molecules below 100 m/s. Again, transverse overfo- 
cusing and longitudinal filtering of the molecular packet 
by the deceleration electrodes are responsible for large 
losses in the decelerator, and this suggests that measures 
beyond modified switching schemes are required to over- 
come these loss mechanisms. 
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FIG. 10: Monte Carlo results of decelerated molecule number 
using S = 3+ and = 20° versus final velocity. The number 
next to each data point is the number of stages used. Because 
of transverse overfocusing, essentially no molecules survive 
below 100 m/s. 
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FIG. 11: (color online) Quadrupole-guiding decelerator. (a) 
Schematic of QGD. (b) Electric field of quadrupole guiding 
stage energized to ±12.5 kV. (c) Switching scheme for decel- 
eration with the QGD. 



VI. QUADRUPOLE-GUIDING DECELERATOR 

In addition to modifying the timing sequences of ex- 
isting decelerators, it is possible, and perhaps preferable, 
to uncouple the longitudinal motion inside the Stark de- 
celerator from the transverse motion by redesigning the 
decelerator electrode geometry. One simple redesign. 
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which we cah the quadrupole-guiding decelerator (QGD), 
is shown in Fig. 11(a). In this decelerator, a quadrupole- 
guiding stage (Fig. 11(b)) is interleaved between each 
deceleration stage. While it may not be necessary to 
have a quadrupole-guiding stage between each deceler- 
ation stage (especially in the beginning of the deceler- 
ator), it simplifies the analysis and will be used here. 
The switching of the electric fields inside a QGD is sim- 
ilar to a traditional decelerator operated with = 1. 
Figure 11(c) shows the potential energy experienced by 
a molecule decelerated at (j)o ~ 45° and is represented 
by a thick black curve. In this panel, the quadrupole- 
guiding electrodes are omitted for clarity. Note that the 
quadrupoles are always energized and their center coin- 
cides with the (j)o = 0° position, while (j)o = 90° occurs 
between the deceleration electrodes. 

Deceleration with a QGD enjoys the same longitudinal 
phase-stability as a traditional decelerator. In a QGD, 
the maximum stable excursion position ^(p^nax 

and ve- 
locity Avmax will differ from that of the traditional de- 
celerator because the decelerating electrodes are most 
likely further apart. In other words, to prevent high- 
voltage breakdown, the quadrupole stages require the 
same inter-stage spacing as deceleration stages in a tradi- 
tional decelerator. Thus, the decelerating electrodes for 
a QGD will be twice as far apart as in our traditional 
decelerator, which possesses an inter- stage spacing of 5.5 
mm. Since the dependence of the decelerating force on 
0o is less steep, the shape of the stable longitudinal phase 
space will change. Understanding the shape of the sta- 
ble longitudinal phase area is crucial for predicting the 
performance of the QGD, and can be derived by exam- 
ining the longitudinal forces inside the QGD. Shown in 
Fig. 12(a) is the on- axis Stark shift of an OH molecule in 
the |2,zb2, — ) state inside the unit cell, defined as 3 de- 



celeration stages of the QGD. The solid line is the Stark 
shift due to the slowing stage centered at 11 mm, while 
the dashed line is the Stark shift of the stages which will 
be energized when the fields are switched. The subtrac- 
tion of these two curves, shown in Fig. 12(b) as a solid 
line, is the amount of energy removed each time the fields 
are switched, AKE. We represent AKE as sum of sine- 
functions [27] 

AKE{(j))= ^ ansin(n0), (1) 

n—odd 

where we have used the definition of the phase angle (j) = 
(2:/L) 180°. A fit of the first three terms of this equation 
to the actual AKE for deceleration stages spaced by 11 
mm is shown as a dashed line in Fig. 12, resulting in 
the fit values ai = 1.221 cm~^, = 0.450 cm~^, and 
as = 0.089 cm~^. Using this fit, we derive the equation of 
motion of the molecules about the synchronous molecule 
position as 



^ + ^(Ai^i?(A</. + <^o)-A7^i?(0o)) = O, (2) 
where we have used the excursion of the molecule from 
the synchronous molecule Acj) = (p — (po- The maximum 
stable forward excursion of a non- synchronous molecule 
is exactly the same as a traditional decelerator and is 
given as [17] 

A<P+,M = 180° -2<Po. (3) 

We calculate the work done in bringing a molecule start- 
ing at this position with zero velocity to the synchronous 
molecule position as 



End rO 

Fdx 



^ (a^[sin(n(A0 + (j)o)) - sin(n0o)])c^(A0). 



I start -^^^taAMn^odds 

Integrating this equation and setting it equal to the kinetic energy yields the maximum stable excursion velocity: 



(4) 



AVmax{(po) 



^ UITT 



where (po is now in radians. 



COS(n0o) .TT , ^ . ^< ^^ 

(2 -0o)sm(0o) ), 



(5) 



r 



Using Eqs. 2- 5 it is possible to solve for the longitudi- 
nal separatrix, which separates stable deceleration from 
unstable motion inside the decelerator. These separatri- 
ces are shown (thick black lines) along with the results 
of Monte Carlo simulations of a QGD in the left column 
of Fig. 13 for successive phase angles. The longitudinal 
phase space is shown with each dot representing the po- 
sition of a stable molecule. The lack of structure inside 



these separatrices is evidence of the lack of coupling be- 
tween the transverse and longitudinal modes. The right 
column, which shows simulated ToF curves, reveals a sin- 
gle stable peak arriving at later times as (j)o is increased. 
These simulations are for a single |2, ±2, — ) state of OH 
and do not exhibit the large background contribution of 
the other states of OH present in experimental ToF data. 
From the comparison of the simulations with the ana- 
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FIG. 12: (a) The Stark shift of an OH molecule in the |2 ± 
2, — ) state inside the QGD. The solid curve is the Stark shift 
due to the slowing electrodes, while the dashed curve is the 
Stark shift due to the electrodes that will be energized at 
the switching time, (b) The change in the molecule's kinetic 
energy as a function of position is shown (solid) as well as a 
fit of Eq. 1, including up to n = 3 (dashed). The solid curve 
is calculated from the subtraction of the two curves in panel 
(a). 



lytical results represented by the separatrices, we see that 
the simple theory of Eqs. 2- 5 is quite accurate in describ- 
ing the longitudinal performance of the QGD. Thus, by 
numerically integrating the area inside these separatri- 
ces, we can predict the longitudinal performance of the 
QGD relative to the traditional decelerator. As seen in 
Fig. 14(a), the energy removed per stage of the QGD is 
less steep with because the decelerator stage spacing 
(in this simulation) is twice that of a traditional deceler- 
ator. For this reason, QGD deceleration with the same 
0o as in a traditional decelerator leads to a faster beam, 
given the same number of deceleration stages. When 
comparing the longitudinal acceptance of the two types of 
decelerators it is important to take the limit of high phase 
angles where both values of AKE converge. Nonetheless, 
the QGD shows significant gain over traditional deceler- 
ation as shown in Fig. 14(b). This gain is primarily due 
to the increased physical size of the stable longitudinal 
phase space resulting from the larger deceleration stage 
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FIG. 13: The left column is stable phase space of molecules 
decelerated inside the QGD. The solid line is the separatrix 
predicted by the theory, while the points represent positions 
of molecule in the 3D Monte Carlo simulations. The right col- 
umn shows the ToF spectra of OH molecules in the |2 di 2, — ) 
state at the exit of this decelerator which has 142 deceleration 
stages, along with 142 quadrupole stages. 



spacmg. 

It is important to note that since Fig. 14 compares 
only the total area inside the separatrix-and S = 1 de- 
celeration does not completely fill this area due to the 
coupling effects-this gain is actually an underestimate of 
the QGD longitudinal performance. Furthermore, these 
graphs do not include transverse focusing effects, which 
can only be properly included through detailed simula- 
tion. The results of Monte Carlo simulations including 
these transverse effects are shown in Fig. 15. The num- 
ber of decelerated molecules versus final speed is plotted 
for both the traditional decelerator operating at 5 = 1 
and the QGD operating at two different quadrupole rod 
voltages, ±1 kV and ±3 kV. While the QGD initially 
delivers more molecules, once the molecules are deceler- 
ated below 100 m/s, the decelerated molecule number 
falls off abruptly. This behavior is expected since, for 
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FIG. 14: (color online) (a) The energy removed per stage as a 
function of phase angle for both traditional deceleration and 
deceleration with a QGD. Both curves are calculated for OH 
in the |2,zb2, — ) state, (b) The calculated phase-stable area 
for deceleration versus phase angle for traditional deceleration 
and deceleration with a QGD is plotted on the left axis. The 
gain of the QGD over traditional deceleration is plotted on 
the right axis. 



these simulations, the voltage on the quadrupole-guiding 
stages is held constant throughout the slowing sequence. 
As detailed in Eq. (3) of Ref. [16], the focal length of a 
transverse guiding element is directly proportional to the 
molecular kinetic energy. Therefore, as the mean speed 
of the packet is decreased, the molecules are overfocused 
and collide with the decelerator electrodes. This can be 
prevented by lowering the voltage on the quadrupole- 
guiding stages during the deceleration process. Figure 16 
displays simulation results of deceleration with this dy- 
namically scaled voltage compared to 5* = 1 slowing at 
(t)o = 30.43°. For this simulation, the quadrupole volt- 
ages are scaled by {v /v initial)^ '^^^ after each decelera- 
tion stage, where v is the instantaneous packet veloc- 
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FIG. 15: (color online) Simulations of the phase-stable 
molecule number as a function of stage number in the QGD 
and S — 1 decelerator. The voltage on the quadrupole stages 
of the QGD is held constant throughout the deceleration se- 
quence. All simulation data is for v final 14 m/s. The 
traces shown are S — 1 deceleration at (j)o — 30.43°, QGD 
operated with ±1 kV on the quadrupoles, and QGD operated 
with ±3 kV on the quadrupoles. Note the decrease in stable 
molecule number in the last several stages for the QGD re- 
sults. This decrease is due to transverse overfocusing of the 
slow molecules through the final few quadrupole stages, and 
suggests that scaling of quadrupole voltage is necessary. 



ity directly following each stage switch. The exponent 
of 0.875 is found empirically to produce the most gain 
at V final = 14 m/s. For ease of simulation, the trans- 
verse forces are scaled by {v / v initial)^ '^^^ whenever the 
molecules are closer to a quadrupole-guiding stage than 
to a deceleration stage. While this may be a poor ap- 
proximation at the lowest speeds, it will likely lead to an 
underestimate of the decelerated molecule number since 
the transverse guidance of the quadrupole-guiding stage 
extends beyond this regime. Even if it leads to an over- 
estimate, proper control of the quadrupole voltages may 
compensate any overfocusing introduced by the deceler- 
ating elements. As seen in Fig. 16, dynamically control- 
ling the voltage of the quadrupole-guiding stages leads 
to a factor of 5 increase in decelerated number for larger 
velocities (>80 m/s) and delivers about 40% more de- 
celerated molecules than traditional S = 1 deceleration 
provides at the lowest final speeds (14 m/s). Because 
the voltages applied to the quadrupole-guiding stages are 
relatively low, dynamic control of them should be possi- 
ble using an analog waveform generator and high- volt age 
amplifier. It should be noted that the optimal voltage 
scaling may vary among decelerators since the real focal 
length depends sensitively on the electrode construction, 
and at low speeds the transverse focusing of the decel- 
erator electrodes becomes significant. This is, presum- 
ably, because the overfocusing of the decelerator elec- 
trodes can be compensated to a certain degree by inject- 
ing molecules which are already slightly overfocused. In 
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FIG. 16: (color online) Monte Carlo simulation results for the 
decelerated molecule number using traditional S = 1 deceler- 
ation (00 = 30.43°) and deceleration using a QGD (00 = 
52.75°) with a dynamic voltage scaling of {v /v initial)^ '^^^ ■ 
For both curves, 142 stages of deceleration were used, and 
V final = 14 m/s. The different phase angles chosen for the 
two decelerators are a result of their differing potential pro- 
files for deceleration. The vertical dashed lines represent the 
deceleration stage at the given velocity. Note that when the 
quadrupole voltage within the QGD is scaled in this manner, 
we observe a 40% gain in molecule number at 14 m/s, and a 
factor of 5 gain over 5* = 1 at higher velocities. 



other words, two focusing stages can overcome the over- 
focusing of a single stage. Thus, it may be possible to 
use adaptive algorithms to optimize the quadrupole volt- 
age or change the design of the decelerating electrodes 
so that they provide less transverse focusing, and maxi- 
mize the number of decelerated molecules beyond what 
is reported here [30]. 

One important advantage of the QGD over traditional 
decelerators is its inherent ability to support more de- 
celeration stages. Because the molecules experience a 
tunable transverse focusing element after each decelera- 
tion stage, there is little loss in efficiency by extending 
the number of deceleration stages. In fact, as seen in 
Fig. 17, there is essentially no loss until the molecules 
are decelerated to the lowest speeds previously discussed. 
This low- velocity loss is due to the aforementioned trans- 
verse overfocusing and longitudinal reflection. The for- 
mer loss mechanism may be overcome, while the latter 
presents a fundamental limit. Even with this loss, the 
QGD outperforms both 6' = 3 (Fig. 8) and 5 = 3+ 
(Fig. 10). Thus, the QGD is an ideal decelerator for more 
efficiently producing cold molecules and, perhaps more 
importantly, the ideal decelerator for slowing molecules 
with poor Stark shift to mass ratios, like H2O and SO2. 



FIG. 17: (color online) 5* = 1 simulation results for v final — 
14 m/s (00 = 5.22°) plotted along with simulation results 
using the voltage-scaled QGD to the same v final {4>o = 23.5°). 
The number labeling each vertical dashed line is the number 
of deceleration stages necessary to reach the given velocity. 
Note the large number of stages (803) used to reach 14 m/s, 
which suggests that a very long QGD may be employed for 
slowing molecules with a poor Stark shift to mass ratio. 



initial experiments and detailed Monte Carlo simulations 
to address them. While the use of decelerator overtones 
yields improvement over = 1 deceleration at high to 
intermediate speeds {v final > 80 m/s), the loss at very 
low velocities remains problematic. The QGD solves the 
problem of coupling between the transverse and longitu- 
dinal motions inside a Stark decelerator by introducing 
independent transverse focusing elements. By dynami- 
cally controlling the focal length (voltage) of these guid- 
ing elements, large improvements (factor of 5) in deceler- 
ation efficiency can be achieved for v final > 80 m/s, and 
gain of ~ 40% is predicted for the lowest velocities. Fur- 
thermore, it appears that with dynamic control of the 
guiding stage focal length there should be no limit to 
the length of decelerator that can be built. This enables 
the deceleration of molecules with a poor Stark shift to 
mass ratio. However, we note that none of the tech- 
niques described in this Article mitigate the longitudinal 
low-velocity loss due to reflection, which appears to be a 
fundamental component of Stark deceleration. Building 
upon the strong correlation between simulation and ex- 
perimental results, we are confident that the simulations 
presented in this work provide a solid guideline for future 
implementations of Stark deceleration. 
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